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Security challenges over the years has led to the need for an improvement in 
the traditional security approaches. This led to the advent of biometrics. 
Recently, among the biometric approaches, sclera has been an area of 
imense study. This is due to its accuracy; however, segmentation of the 
sclera has been a limiting factor to the application of this biometric trait. 
Several approaches have been proposed in literature but there is still the need 
to improve the segmentation accuracy. This study proposes the use of 
circular hough transform and a modified run-data based algorithm. The 
study also presented a sclera recognition system using the compound local 
binary pattern for features extraction and Manhattan distance for 
classification. The system produced a segmentation accuracy of 99.9% for 
sclera blood vessels, periocular and iris (SBVPI) sclera database and 100% 
for manually captured sclera database. The system produced an accuracy of 
99.98 for SBVPI sclera database and 99.99% for manually captured sclera 


database. 
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1. INTRODUCTION 

The increase in the growth of digital technology has led to the need for better security methods. This 
led to the increase in the application of biometric technology in identity verification and authentication. This 
increase is due to the reliability and unique nature of biometric traits. Biometric authentication uses a 
person’s behavioural or physical trait for recognition [1], [2]. These traits include iris, fingerprint, palmprint, 
iris, signature, and among other [2], [3]. These traits can be categorized as being physical or behavioural. The 
physical trait answers the question of who an individual is, while the behavioural trait answers how an 
individual does something. Among the physical trait are those that are related to an individual’s ocular 
properties (that is the iris, pupil, and sclera) [1], [3]. 

The iris is noted for its accuracy in near infrared images; however, this drops in visible light 
wavelength with no constraint during capturing [4]. To alleviate the challenges associated with both traits, 
pairing the trait with another would improve the performance of iris [5], [6]. This is termed multimodal 
biometrics; the use of more than one biometric trait for recognition. This is usually aimed at increasing 
accuracy and increasing the difficulty of tricking the system. Among the other suitable biometric traits that 
can be combined with iris, sclera seems to be most suitable [6], [7]. This is because it can be captured with 
the iris. Iris is the outer circle of the human eye, while the sclera is the white part of the eye. 
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Sclera trait is an area of biometric that has recently been receiving attention and this is due to the 
fact that it is a highly-protected part of the eye that is difficult to forge. The vessel pattern on the sclera is 
studied and found to not be the same for two people and different for two pairs of eyes [8], [9]. This 
uniqueness makes it suitable for use as a biometric trait and there is a need to develop an efficient recognition 
system for the trait. However, the segmentation of the sclera from the eye image has limited the application 
of this biometric trait. Several studies have proposed different approaches for sclera segmentation; however, 
the segmentation accuracy shows that there is still a lot of room for improvement [10], [11]. Hence, this 
study examines the segmentation of sclera from the eye image. It also presents a sclera recognition system. 


2. RELATED WORKS 

Sclera trait in biometric is important in biometric as it can uniquely identify an individual. It can 
further be combined with other ocular biometrics to increase the accuracy of the system. Recent trends have 
examined the segmentation of sclera in the eye image and use of this biometric trait for recognition in a 
unimodal and multimodal way. Literatures relevant to both sclera segmentation and recognition in this 
section. 


2.1. Sclera segmentation 

Rot et al. [1] developed a cascaded convolutional neural network (CNN) assembly for segmenting 
sclera. They went further to propose a CNN model (ScleraNet) for extraction of discriminatory features from 
the sclera. A new sclera dataset of ocular images named sclera blood vessels, periocular and iris (SBVPI) was 
collected. An accuracy of 93.3% was recorded. In the study by Xu et al [2], they proposed a model for sclera 
segmentation. It starts with filtering using median filter. The Cr colour channel image. Saviola method for 
iris segmentation. Global otsu algorithm and adaptive otsu algorithm were applied to the two segmented 
images and the output was fused together. the local area histogram equalization was used for image 
enhancement. The length of blood vessels was used for discrimination. They recorded a segmentation 
accuracy of 98.1%. 

Pathak et al. [3] also proposed sclera segmentation system. Their system starts with preprocessing, 
using bilateral filtering. After this, features were extracted in the form of brightness, colour and texture. 
Using the extracted features, CNN was used for the segmentation of the sclera, iris and pupil based on the 
entropy value. An accuracy of 98% and 99.423 was obtained on Multimedia University (MMU) and 
Universidade de Beira Iris (UBIRIS) v.2 datasets respectively. 

Maheshan et al. [4] proposed the segmentation of sclera from the eye image using a modified 
intuitionistic fuzzy clustering approach. This is aimed at improving the performance of the traditional fuzzy 
set that non-membership value is always the complement of the membership value. They obtained a precision 
of 85.10 and a recall of 77.90. In the study by Vitek et al. [5], sclera segmentation was performed using a 
region growing method called Active contour without edge. In their approach, two points are initialized, after 
which the active contour without edge is performed. 


2.2. Sclera recognition 

Among the recent study in sclera recognition is the work of Pathak et al. [6]. The paper presented a 
multimodal biometrics recognition using the pupil, iris and sclera. The preprocessing starts with the use of 
min-max normalization. After which entropy based CNN was used to segment the image. The colour 
histogram features, gabor features and Y-shaped features were used to calculate the support value of each trait 
respectively. UBIRIS V2 database was used for testing the system. An accuracy of 97.99% was recorded. 

In the study by Maheshian et al. [7], they explored the use of deep learning for sclera recognition. In 
particular they proposed a custom layer CNN that they termed CNN sclera recognition engine (CNNSRE). 
The network model is made up of four convolutional units and a fully connected layer. An accuracy of 87.65 
was recorded. 

Lee and Kim [8] proposed a sclera recognition system using bifurcation-based descriptor. The sclera 
was segmented from the ocular image, and the blood vessels were extracted. Morphological operation was 
used to thin the blood vessels extracted. The bifurcation of the blood vessels was extracted and the distance 
between the central bifurcation and its neighbours were used as features. Each local structure is then 
matched. The result showed an equal error rate (EER) of 1.92%. 

The study by Zhu et al. [9] proposed a sclera recognition method that consist of a non-learning 
based segmentation method and a supervised learning technique for classification. Contrast limited adaptive 
histogram equalization (CLAHE) was applied to segment the sclera from the eye image. Their approach 
introduced a stem-and-leaf branches network for classification. A new sclera dataset named ScleraVO was 
introduced. Their system showed an accuracy of 96.91. 
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3. METHOD 

The block diagram of the system is shown in the Figure 1. This details the key operations to be 
performed by the system. It includes collection of eye images, filtering of each eye image, segmentation, 
features extraction, and matching. 
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i | 
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Figure 1. Block diagram of the designed system 


3.1. Eye image 

For this system, the SBVPI eye database and manually captured eye images were used [1]. For the 
SBVPI database, 1000 eye images were used for testing. 500 eye images were genuine and the remaining 
half not genuine (not registered on the system). For the manually capture eye images, a canon EOS 60D 
camera was used to capture 200 eye images. Both eye image dataset was used in the developed system. A 
sample SBVPI eye image is shown in Figure 2. 


Figure 2. A sample eye image 


3.2. Grayscale conversion 

This involves representing the image in its gradient form. Computing the average of the red, green, 
and blue value for each pixel gives the gradient form of the image. Figure 3 shows the grey output of 
Figure 2. 


Ne ae 
a © Š aine Fate 


Figure 3. Gray eye image output 
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3.3. Canny edge detection 
Canny edge detection algorithm is used to convert the image to a binary form which depicts the 
edges of the image. It consists of the following steps [10], [11]. 


3.3.1. Filtering 
Filtering is aimed at reducing noise and unwanted textures in the ocular image. In (1) and (2) shows 
the image f (m,n) filtered into g(m, n). 


g(m,n) = G_o (m,n) * f(m,n) (1) 
where 
G_o (m,n) = 1/V( [no] *2) exp(—(m*2 + n*2)/ [20]^2) (2) 


f (m,n) represents the image. 
Where m is the distance in the horizontal axis, n is the distance in the vertical axis, and o is the 
standard deviation of the Gaussian distribution [12]. 


3.3.2. Finding the gradient 

The gradient measure on the image is performed with the sobel operator. The edge detection 
operator (sobel) returns a value for the first derivative in the horizontal direction (gm) and the vertical 
direction (gn). From this the edge gradient and direction can be obtained using (3) [13]. 


M(m,n) = V( lgm} 2 (m,n) + [gn] 2 (m,n)) (3) 
The sobel operator estimates gradient in the m and n direction using two convolutional masks. The 


masks are shown in Figure 4. The gm estimates gradient in the m direction and the gp estimates gradient in 
the n direction. 


Figure 4. Two tables showing sobel masks (kernel), gm and gn respectively 


Where the direction of the edge was computed with (4) using the gradient in the m and n direction, 
A(m,n) = [tan] *(—1) [g_n (m,n)/g_m (m,n)] (4) 


3.3.3. Edge direction 

After determining the edge direction, the obtained edge direction is then converted to a direction that 
can be traced in the image [14]. The possible direction that can be traced can one of the eight neighbouring 
pixels. The eight directions are north, north-east, north-west, right, left, south, south-east and south-west of 
the pixel being considered. 


3.3.4. Suppress non-maxima pixels 

After the direction of the edges are gotten, non-maximal suppression is performed on the image. To 
trace along the edges and set to zero pixels not considered edges, non-maximal suppression was used. The 
output of this step is a thinned edge line [15]. 


3.3.5. Double thresholding 

Two thresholds were chosen as T1 and T2 (that is T2 is the major edge and T1 is a minor edge) such 
that (T1<T2) to obtain the binary image. Edge segments of T2 were linked to obtain the continuous edges 
and any that was below T2 and was greater than T1 that lies along the continuous T2 edge segments was 
consider a major edge [16]. 
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3.3.6. Linking of edges 

The final step involves linking the edge segments to form continuous edges [11], [14], [15]. Edge 
image from canny edge detection algorithm as shown in Figure 5. Figure 5 shows the edge image of 
Figures 2 and 3. 


Figure 5. Edge image from canny edge detection algorithm 


3.4. Circular hough transform 

Circular hough transform is used to detect the iris boundary and the pupil boundary. The output of 
this gives the center coordinates of the iris and the pupil [3], [17], [18]. It also gives the radius of the iris and 
the pupil. The CHT is represented in (5) [12], [19]. Figure 6 shows the pupil boundary detected after using 
the hough transform on Figure 5. 


((x-a)] 24+ Ky- b)J ^2=r^2 (5) 


With a and b being the center of the circle in x and y directions respectively and r being the 
radius. The parameter of the circle is as shown in (6) and (7) [20]-[22]. 


x =a + rcos(8) (6) 


y =b +rsin(0) (7) 


Figure 6. Pupil boundary detection using circular hough transform 


3.5. Segmentation using the modified run-data-based following algorithm 

The traditional edge-point tracing of run data is used in the run-data-based following method, which 
employs run data in pairs consisting of an object's left and right edges that are produced using horizontal scan 
lines from left to right on an image. There may be an outer contour and multiple inside contours on the 
object. There are five types of run data: left edge of outer contour, right edge of outer contour; left edge of 
outer contour, left edge of inner contour; right edge of inner contour, left edge of inner contour; right edge of 
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inner contour, left edge of inner contour; right edge of inner contour, right edge of outer contour; right edge 
of inner contour, right edge of outer contour; and right edge of outer contour, left edge of outer contour [23]. 

However, for extraction of sclera, modifications were proposed to the traditional run-data-based 
algorithm. The modified run-data-based following algorithm was applied after the detection of the iris 
boundary. The key-points were located using the center co-ordinates of the iris boundary and the radius of the 
iris that were detected. The key-points detected are the eight points at the right sclera and the left sclera. 
These points helped in cropping out two images from the left and right sclera. To locate the points, a 
modified run-data-based following (MRDBF) algorithm was used. It must be noted that the MRDBF 
algorithm does not directly trace the boundary but picks points (8 points) that forms two rectangular 
boundaries within the sclera regions (right and left) as shown in Figure 7. The changes made to the algorithm 
include: the starting point for the algorithm is dependent on the circular hough transform. Also, the tracing is 
done within the boundary of the image with the purpose of locating distinct and consistent points within the 
boundary and this is helpful in the case where the boundary is not completely defined. As opposed to 
scanning the image row-wise and column-wise to detect the inner and outer pixels, the algorithm scans the 
image from the reference points. The MRDBF algorithm is shown in Algorithm 1. 

Figure 8 shows the flow chart of the proposed algorithm. The points detected in Figure 7 are shown 
in Figure 9. Using the eight key-points (P1-P8), the sclera is cropped out. Figure 10(a) shows the sclera 
image cropped out using eight key-points on the left side of the iris. Figure 10(b) shows the sclera image 
cropped out using eight key-points on the right side of the iris. 


y-axis 


Note: Where x & column and y is row 


cor dinate 


0 X-axis 


Figure 7. Modified run-data-based following algorithm 


Algorithm 1: Modified Run-Data-Based Following Algorithm 
Input: An array of binary pixels of an eye image (P) with P={0,1}. Where 0 is the 
background pixel and 1 is an edge pixel. 
Let mandn be the dimension of the input array respectively. 
Let CHT_row and CHT_col be the row and column center co-ordinate of the iris returned by 
the CHT (Circular Hough transform) 
Let Iris_r be the radius of the iris. 
Start: 


1. xplotl = CHT col + Iris _r + 4 
2. yplotl = CHT_row 
3. for b = xplotl:1: n (i.e. increment by 1 and stop at end of column) 
if P (yplotl, b) == 1 (i.e. checking for foreground) 
xplot2= b (i.e. detect the column point for xplot2 (length of the 
sclera region) ) 
break 
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else 
xplot2 = n 
end if 
end 
4. xplot3 = CHT col ~ Iris r - 4 
5. yplot3 = CHT_row 
6. for b = xplot3: -1: 1 (i.e decrementing by 1) 
if P (yplot3, b) == 1 (i.e. checking for foreground) 
xplot4 =b 
break 
else 
xplot4 == 
end if 
end 


//Detecting points for the first rectangle 
7. xplot_new_poin = (xplot2- xplotl)/ 2 
8. xplot_new point 

and xplot 2 for a rectangle) 
9. yplot_new_point = yplotl 


(xplotl + xplot_new_point) (i.e. getting a point between xplotl 


10. for y=yplot_new point: 1: m (i.e. incrementing by 1 to the end of row) 


if P (y, xplot_new_point) == 1 
y2plot = y (i.e. new point y2plot) 
break 
end if 
end 
11. for y = yplot_new_ point: -1: 1 (i.e. decrement by -1) 
if P (y, xplot_new_point) == 1 
yplot = y 
break 
else 
end if 
end 


//points detected are: 

//Pl(yplot, xplot_new_point) 

//P2(y2plot, xplot_new_ point) 

//and 

//P3(yplot, xplotl) 

//PA(y2plot, xplotl) 

// Detecting points for the left rectangle 
12. xplot_new_poinL = (xplot3 - xplot4)/2 
13. xplot_new_pointL = (xplot3-xplot_new_point) 
14. yplot_new_pointL = yplot3 
15. for y = yplot_new_pointL: 1: m 


if P (ya xplot_new_pointL) == 1 
y2Lplot = y (i.e. new point y2Lplot) 
break 

else 

endif 


end 


16. for y = yplot_new_point: -1: 1 


if P (y, xplot_new_pointL) == 1 
yLplot = y (i.e. new point yLplot) 
break 

else 

end if 


end 


//points detected are: 
//P5(yLplot, xplot_new_pointL) 
//P6(y2Lplot, xplot_new_pointL) 
//and 

//P7(yLplot, xplot3) 
//P8(y2Lplot, xplot3) 


Output: The points detected are the points Pl, P2, P3, P4, P5, P6, P7, 
respective co-ordinates. 
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Edge 
Image (P) 


Detect point Xplot1 on the 
horizontal axis (i.e., Xplot1 = 
CHT _col + Iris_r+4), from P 


Detect point Xplot3 on the 
horizontal axis (i.e., Xplot3 = 
CHT_col - Iris_r+4), from P 


Scan edge image P to detect 
Xplot2 using incrementing 
loop 


Scan edge image P to detect 
Xplot4 using decrementing 


Scan edge image P to detect 
point by averaging Xplot1 


loop and Xplot2 
Scan P vertically from Scan P vertically from 
Xplot_new_point to detect Xplot_new_pointL to detect 
yplot and y2plot. yLplot and y2Lplot. 


Detect Points: P1(yplot, xplot_new_point) 
P2(y2plot, xplot_new_point) 
P3(yplot, xplot1) 
P4(y2plot, xplotl) 
P5(yLplot, xplot_new_pointL) 
P6(y2Lplot, xplot_new_pomtL) 
P7(yLplot, xplot3) 
P8(y2Lplot, xplot3) 


Combine coordinate to locate 
points above and below 
Xplot3 and Xplotl 


8 points 
detected 


Crop out two sections from 
edge image P with the 8 
points detected: (P1,P2,P3 
and P4), (P5,P6,P7 and P8) 


Not an Eye image 


Figure 8. Flowchart of the modified run-data-based following algorithm 


Figure 9. Figure showing the eight key-points (P1-P8) 
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(a) (b) 


Figure 10. Sclera image extracted from the of the eye image; (a) sclera image of the left side of the iris and 
(b) sclera image of the left side of the iris 


3.6. Features extraction 

The feature extracted from the sclera image is the gradient of the segmented sclera regions. In terms 
of the segmented iris image, the gradient (patterns) of the iris was extracted. For the features extraction of the 
iris and the sclera, compound local binary pattern (CLBP) was used. LBP is known for its computational 
simplicity and resistance to changes in illumination; however, it discards the magnitude information of the 
difference between the center and the neighbour gray values in a local neighbourhood. This may result in 
inconsistent codes which is resolved by CLBP. The equation for CLBP is shown in (8) [24], [25]: 


fe in —ig <0, lip — ic] < Mang 
01 t,-i, <0, |ip— icl > Mang 
oe in —te 20, [ly — tc] S Mavg 

11 Otherwise 


S(ip, ic) = 


(8) 


where, i, is the grey value of the center pixel, i, is the gray value of a neighbour p, and Mayg is the average 
magnitude of the difference between i, and i, in the local neighbourhood. Figure 11 shows the equivalent 
CLBP of the extracted sclera images after they were joined together. 


Figure 11. CLBP of sclera image 


3.7. Matching 

For the matching, Manhattan distance was used. Manhattan distance is the sum of the absolute 
difference between two vectors in n-dimensional vector space. According to Shang et al. [26], it is the sum of 
the absolute differences of their corresponding components. It is as shown in (9) [27]. 


d(X,Y) = Xi- lxi — vil (9) 


3.8. Performance evaluation metrics 

The evaluation of the developed system was performed using the false acceptance ratio (FAR), the 
false rejectance ratio (FRR), and the accuracy. The methods were selected because of the frequent use in 
literature. The FAR is the ratio of falsely accepted users to the total false users tested. The FRR is the ratio of 
falsely rejected users to the total genuine users submitted. The FAR and FRR are also known as the false 
positive ratio (FPR) and false negative ratio (FNR) respectively. The accuracy of the system is the percentage 
of the users that were correctly identified by the biometric system. To compute the FAR and FRR, the true 
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positive (TP), false negative (FN), false positive (FP), and true negative (TN) are required. The TP is the 
number of genuine individuals that were positively predicted. The FN is the number of genuine users that are 
predicted negative. The FP is the number of impostors that were predicted positive. The TN is the number of 
impostors that were negatively predicted. In (10) to (12) shows the FAR, FRR and accuracy [28], [29]. 


FAR = FPR = FP/(FP + TN) (10) 


FRR = FNR = FN/(TP + FN) (11) 


Accuracy = 100 — “EHR (12) 


4. RESULTS AND DISCUSSION 

Tables 1 and 2 show the confusion matrix of the manually captured eye images and the SBVPI eye 
image database. From Tables 1 and 2, (10) and (11), the FAR and FRR obtained after testing the system are 
both 0.01. This produced an accuracy of 99.99 for the manually captured eye images using the canon camera. 
The SBVPI database of eye images gave an accuracy of 99.98, with the FAR and FRR being 0.016 and 0.024 
respectively. The EER for the system was used to select the optimum threshold for matching. The EER for 
this system is shown in Figure 12. 


Table 1. Confusion matrix of the manually captured Table 2. Confusion matrix of SBVPI eye image 
eye images database 
Total no=200 Predicted-no _ Predicted—yes Total no=1000 Predicted-no _Predicted—yes 
Actual—no TN=99 FP=1 100 Actual—no TN=488 FP=12 500 
Actual-yes FN=1 TP=99 100 Actual-yes FN=8 TP=492 500 
Total 100 100 200 Total 496 504 1000 


EQUAL ERROR RATE (EER) 


FAR 


THRESHOLD 
BR 


FAR / FRR 


Figure 12. EER of the system 


4.1. Comparison of performance 

The comparison of the accuracy of the segmentation algorithm is presented in Table 3. The 
segmentation algorithm proposed in this document achieved an accuracy of 99.9% using the SBVPI eye 
image database. The segmentation accuracy recorded for the manually captured eye images was 100%. This 
accuracy shows a higher performance without the required training for deep learning techniques when 
compared with the closest method in Table 3. A close look at the comparison of the system’s segmentation 
performance is shown in Table 4. 


Table 3. Segmentation accuracy comparison table 


Paper Method Accuracy (%) 
Rot et al. [1] Saviola method SBVPI 98.1 
Dong et al. [30] Otsu thresholding UBIRIS v1 98.1 
Pathak et al. [3] CNN MMU 98 
UBIRIS v2 99.423 
This study Circular hough and modified run-data based following algorithm SBVPI 99.9 


Manual capture 100 
Table 4. Comparison of the proposed system performance with other systems 
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Paper Accuracy (%) 
Pathak et al. [3] 97.99 
Maheshan et al. [7] 87.65 
Proposed 99.98 (SBVPI) 


99.99 (Captured database) 


5. CONCLUSION 

The advent of biometric traits for solving security challenges has been an area of much research. 
More traits are being discovered, amongst which is the sclera. However, the segmentation of the sclera in an 
eye image has been an area that requires more accuracy. This study has proposed an approach for the 
segmentation of the sclera. The approach uses circular hough transform and a modified run-data based 
algorithm. The accuracy of the presented segmentation algorithm showed an improvement when compared 
with previous approaches. This sclera segmentation technique is robust as it allows for the extraction of 
sclera on both sides of the iris. The proposed algorithm also made use of circular hough transform to 
recognize the iris. This can make the segmentation of both the sclera and the iris seamless. It also requires 
less computation resource and no training when compared to supervised learning-based techniques. The 
paper presented a sclera system using the proposed segmentation technique and compound local binary 
pattern. The system showed an encouraging performance. 
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